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I. INTRODUCTION 



Both the massive lepton pairs production in hadronic colhsions (Drell-Yan) and deep 
inelastic scattering (DIS) at high energy are the most outstanding processes probing the 
hadron structure. On the deep inelastic side, a large amount of work has been done to 
describe the copious data at medium and very small Bjorken scaling variable x, based on 
perturbative QCD In the high energy domain it has been found that important unitarity 
corrections should be taken into account regarding the standard pQCD approach . These 
phenomena are currently denominated perturbative shadowing or saturation effects 
Concerning the Drell-Yan sector, the pQCD tools have produced a reasonable theoretical 
understanding of the main observables, despite the small data set available at present 0. 
The forthcoming accelerator experiments (RHIC and LHC) will scan the high energy limit 
of the hadronic reactions and open a new kinematic window, i.e. smaller x values. In 
particular, the quark-gluon plasma (QGP), a new state of the hadronic matter predicted by 
QCD, is expected to be found there p]]. The theoretical description of the QGP production 
is directly associated to a complete knowledge on saturation effects and the transition region 
to the high parton densities. In a specific way, since the production scheme for J/\E' is similar 
to the Drell-Yan one and the latter does not contain final state effects, DY can be considered 
as a baseline process to study J/\E' suppression as a signature of QGP formation [^. 

In the fast proton system, the QCD factorization theorem leads to describe the hadronic 
processes through the convolution of the parton distribution functions (pdf's) with the 
partonic subprocesses. The latter are completely calculated in pQCD up to higher orders, 
whereas only the evolution in the factorization scale of the pdf's is determined. Namely, 
the parton distributions are solutions of the DGLAP evolution equations, whose formalism 
has been successful in describing both DIS and DY data Recently, an alternative 

way to study electron-proton and hadronic reactions is claimed by the color dipole picture 
considering the rest frame description based on /c^-factorization [Q,|l^,|Tl|] . Thus, the basic 
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blocks are the dipole light-cone wavefunction and the dipole-target cross section. Such 
an approach has produced an unified way to study the mentioned processes, however its 
complete connection with the standard DGLAP formalism is not provided yet and deserves 
further studies. 

In the infinite momentum frame, the DY process corresponds to the annihilation of a 
quark (antiquark) from the projectile with antiquark (quark) of the target into a virtual 



photon (vector boson), which afterwards decays into a lepton pair |jT2|. In the leading order 
(LO) calculation, the DY process has a simple electromagnetic character and it can be 
promptly given by QED theory. However, the perturbative QCD results at higher orders 
modify this simple picture. At present, pQCD calculations have been developed up to 



the second order of the strong coupling constant |]T3[. For practical considerations, in 
general the involved next orders contributions are taken into account by a phenomenological 
parameter, namely a K factor which is dependent on the DY kinematic variables. 

In the rest frame, the DY process looks like a bremsstrahlung of a virtual photon decaying 



into a lepton pair, rather than a parton annihilation |TT|]. The bremsstrahlung of the virtual 
photon can occur after or before the interaction with the gluonic field of the target. The 
advantage of this formalism is that the corresponding cross section can be considered in 
terms of the same dipole cross section extracted from small-x DIS in the color dipole picture 



IJ] . At high energy, the unitarity corrections should be included in the dipole cross section. 
Such effects have been considered, for example, in the phenomenological model of G.Biernat- 
Wusthoff (GBW) ||15[, which describes DIS and ep diffractive process with good agreement. 
We notice, however, that the unitarity corrections to the inclusive observables, i.e. total cross 
section or F2, can be hidden into the parametrization based on DGLAP approach, absorbed 
in the initial conditions, thus providing an excellent data description as seen in the updated 



NLO QCD fits |]T6[. More exclusive observables should be useful to clarify this important 
aspect. After this short remark, we proceed our argumentation. The main disadvantage 
in GBW is that a dynamical explanation of the saturation phenomenon is lacking. On the 
other hand, the Glauber-Mueller approach provides a theoretical development concerning 



parton saturation 0], constraining the pQCD description of the dipole cross section. Here, 
we make use of this formahsm to perform a description of the DY process in the rest frame. 

The goal of this work is to perform a study of DY at high energies considering the color 
dipole picture, in a similar way of recent works [0. Our contribution is based on the use 
of the dipole cross section calculated in perturbative QCD, through the Glauber-Mueller 
formula [Q, which encodes the unitarity effects (saturation) in the parton densities. This 
approach takes into account the multiple Pomeron scattering hypothesis in an eikonal way 
keeping the unitarity of the considered process. A comparison between the phenomenological 
GBW dipole cross section and the theoretical Glauber-Mueller one is presented, verifying 
that the two approaches have different behaviors at higher energies. This is due to the 
dynamical dependence on the gluon distribution in the Glauber-Mueller approach, whose 
Born term recovers the DGLAP kernel in double log approximation (DLA). The nonpertur- 
bative region, i.e. large dipole sizes contributions, is addressed considering the freezing of 
the gluon distribution under the initial perturbative evolution scale Qq. Then, we present 
DY calculations in the rest frame of the target at leading order in a pp collision and perform 
a comparison with the low x DY differential cross section from the E772 Collaboration | |18[| . 
We also produce estimates for the cross section at RHIC energies. 

The outline of this paper is the following. In the next section we present a brief review of 
the DY process in the dipole color picture, discussing the range of validity for this approach 
and showing the role played by the 7*g wavefunction. In Sec. (3) we present high parton 
density effects calculated from the Glauber-Mueller approach, and confront them with the 
phenomenological GBW model. We estimate the contribution of the saturation effects for 
the dipole cross section in high energies (LHC and RHIC). In Sec. (4), a parameter-free 
prediction to the differential DY cross section for the available data at small-a; and estimates 
to RHIC are performed. Finally, in the last section the results are discussed and we present 
our conclusions. 
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II. DRELL-YAN IN THE COLOR DIPOLE PICTURE 



Before tlie description of the Drell-Yan process in the rest frame, we would like to review 
the main kinematical variables and the standard calculations in the laboratory system. 
This is important to clarify the connection between them and to emphasize the asymmetry 
projectile-target in the rest frame picture. 

In the laboratory system, the lepton pairs are produced in the Drell-Yan reaction where 
partons from the projectile (fast proton) interact with the proton target |[T^. Looking at 



the parton level, a quark-antiquark pair annihilates into a virtual photon in leading order 
gg — i> 7* — s> The symmetry between target and projectile is very clear, namely we 

cannot distinguish a quark coming from the proton target or from the incoming beam. The 
momentum fraction carried by the quark from the projectile is labeled Xi and from the target 
is X2- The partonic subprocess above is well known from QED, and the hadroproduction cross 
section is obtained folding the partonic cross section with the quark (antiquark) densities 
evaluated at the invariant M^, the squared lepton pair mass, chosen here as the factorization 
scale yUfa^j,. Their evolution in is given by the standard DGLAP evolution equations. 
Therefore, the DY differential cross section in leading order is given by 
dV'"' 4iraL 1 



where q[q]f{x, M"^) are the corresponding quark (antiquark) densities with flavour / and 
squared charge ej-. The center of mass energy squared is s and the usual notation is 

Xf = Xi- X2, (2) 
T = XiX2 = M'^/s. (3) 

The momentum fractions are rewritten as 

= ^{^xI + 4t + xf), (4) 
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X2 = ^{\/xI + 4t-xf), (5) 



where xp is the longitudinal momentum fraction, labeled Feynman x. Indeed, xp, M and 
s are the kinematic variables experimentally measured, whereas the partonic variables, Xi 
and X2, are reconstructed from them. 

When we consider the target at rest, the DY process looks like a bremsstrahlung: the 
quark from the projectile radiates a photon, which carries a fraction a of the light-cone 
momentum of the initial quark, later decaying into the lepton pair [ see Fig. (1)]. The 
interaction with the target can occur before or after the photon emission. Thus, although 
diagrammatically no dipole to be present, the interference among graphs results in a product 
of two quark amplitudes in the DY cross section, testing the external gluonic field at two 



different transverse positions [0. Therefore, a remarkable feature emerging is that the 
7*g — interaction can be described by the same dipole cross section as in DIS [1^ . 

In the {a, r±) mixed representation, the photoabsorbtion cross section in deep inelastic 
is described by the convolution of the wavef unctions, from the virtual photon and 

the interaction dipole cross section, cr^g. The wavefunctions are considered taking into 
account the first photon Fock state configuration, namely a qq pair. The dipole cross section 
is modeled phenomenologically based on a matching between the hard and soft pieces, 
constrained by the DIS available data. The transverse separation of the qq pair is r±, 
and each quark (or antiquark) of the dipole carries a momentum fraction a (or 1 — a), from 
the incoming photon. The small dipole size configurations can be described through pQCD, 
whereas the large size ones belong to the nonperturbative domain. Hence, one can write the 
photoabsorption cross section as a function of the scaling variable x and photon virtuality 

in the quantum mechanics form [p!0|] . 
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FIGURES 




FIG. 1. The Drell-Yan process in the rest frame, depicting one of the possible interactions 
7*(7-target (see text). 



(6) 



where T, L indicate the transverse and longitudinal contributions to the total cross section. 
In a similar way, the cross section for radiation of a virtual photon from a quark after 
scattering on a proton has the following factorized form in the color dipole picture 

dcrT,L{qp qi*p) 



d\n.a 



(7) 



where we have the same dipole cross section as in DIS. Here r± is the photon-quark transverse 
separation, ar^ is the qq separation and a is the fraction of the light-cone momentum of the 
initial quark taken away by the photon. We notice the difference with the DIS case, where 
the dipole separation is just rj_. Here, aqq is the cross section for scattering a qq pair off a 
proton which depends on the qq transverse separation, and which should take into account 
the saturation effects at high energy. 

The physical interpretation of (|^) is similar to DIS in the light-cone approach (LC). The 
projectile quark state is expanded in its Fock space in the form ||11|], 



|g) = Z2|g) + <^|g7*) + ... 



(8) 



where here one has the expansion in terms of the eigenstates from the quark projectile. 
Instead, in deep inelastic the expansion is constructed from the eigenstates of the incident 
photon [jlO|. Here Z2 is the renormalization constant. 



The well known LC wavefunctions can be calculated in perturbation theory |TyjTT|, and 



depend on the transverse separations and momentum fraction a. They play an important 
role in the dilepton mass M dependence. We take the same notation for the LC wavefunc- 
tions from [ p^ , 



^^.,(«, = ^ (m}a'K',ir^r^) + [1 + (1 - af] v^K^Vr^)) , (9) 



^^,(«, r±)r = - c^fKiivr^) . (10) 



The functions Kq and Ki are the modified Bessel functions and the auxiliary variable r), 
depending on the quark mass m^, is given by 

T]"^ = {l-a)M^ + a^mj. (11) 

The hadronic differential cross section for the Drell-Yan process is expressed in a factor- 
ized form, embedding the partonic cross section, Eq. (1^), into the hadronic environment, in 



the following way [0 



da^^ xi da 



dM'^ dxp GvrM^ (xi + X2) Ai a 



E 



2 



a / \ a 



da{qp — > qj*p) 
din a 



;i2) 



where Cj is the quark charge. In this frame we use standard kinematical variables Xi = 
{2P2.q)/s and X2 = (2Pi.g)/s, with X1X2 = (M^ + qT)/s, where Pi, P2 and q are the four 
momenta of the beam, target and virtual photon, respectively. = and q^ are the 
dilepton invariant mass squared and the squared transverse momentum, respectively. 

The frame dependence of the space-time interpretation of the DY process can be illus- 
trated by different meanings of Xi in different reference frames: we know that in the Breit 
frame, Xi is the momentum fraction of the projectile quark (antiquark) annihilating with 
the target antiquark (quark). In contrast, evaluating the scalar product referred above in 
the target rest frame shows that the projectile quark carries momentum fraction x = Xi/a 
(which is larger than xi), of the parent hadron, and correspondingly, Xi is the momentum 
fraction of the proton carried by the photon. The variable X2 is the momentum fraction of 
the proton carried by the gluon exchange in the t-channel. We are benefited with the fact 
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that the parton densities qi and qi of the projectile enter in the combination , which is 
the structure function of the proton. Therefore, we can rewrite the equation above in the 
following way, 

da^^ 1 da „ f Xi\ dcr^qp q'~i*p) 



dM'^dxp QtiM"^ {xi+X2)Jxi a \aj d\na ' 

where the summation of the longitudinal and transverse contribution was considered. The 
factor aem/(6 vrM^) is due to the photon decay into the lepton pair, coming from electrody- 
namics, the differential cross section da{qp q'y*p)/dlna is taken from equation (|^) and 
our input to aqq in this work is given by the standard gluon distribution in the target 
corrected by saturation effects in the high energy limit. In Eq. (p!3D, the structure of the 
projectile is described by the -F^x, Q^) structure function. 

In the rest frame, the process is asymmetric concerning projectile and target, in contrast 
with the symmetric picture in the Breit frame. The dipole color picture is valid for small X2 
and it takes into account only the gluonic (sea quarks) sector from the target, disregarding 
its valence content. However, both valence and sea quarks in the projectile are parametrized 
in the proton structure function in Eq. (|T^) (for a complete discussion see Refs. [0|). 
Although at present there is little range of experimental measurements in the kinematical 
limit of validity of the color dipole approach, it should provide reasonable results when one 
considers smaller X2 than the currently available. The high energy accelerators LHC and 
RHIC will open a wider kinematical window towards smaller X2 values allowing to test rest 
frame calculations properly. 

To conclude this section, we analyse the behavior of the wavefunctions in the relevant 
kinematic variables. As it will be shown, they play the role of a weight to the dipole cross 
section concerning the transverse separations. In Eqs. (^^, large r± configurations are 
suppressed in the integrated cross section, controlling the nonperturbative contributions 
(large transverse distances domain) to the observables. In the deep inelastic case, the LC 
wavefunctions dependence on the radius at fixed photon virtuality is discussed in 



Refs. ||T^. For the Drell-Yan case, the weight functions are given by: 



Wli{r^,M^) = r±J^ F^{x,/a,M') |v^^4(a, r^)^ • (14) 

In Fig. d^) one shows separately the longitudinal and transverse results for W^'^{rj_, M^) 
as a function of the photon-quark transverse separation r_|_ at fixed lepton pair mass M. The 
choosen momentum fraction was Xi ^ xp = 0.525, since it is a typical experimental value 
(see Sec. 4). Considering this Xi value, the proton structure function is insensitive to the 
lepton pair mass range because it is in the scaling region. Regarding the quark mass, here 
one takes an effective light quark mass = 0.2 GeV in the wavef unctions. 

For the transverse contribution, meaning the upper plot in Fig. (0), one verifies that the 
weight function selects from small up to intermediate photon-quark sizes. This means that 
it is selecting small dipole sizes {ar±) in a similar way to deep inelastic scattering, since 
Xi < a < 1. For our purpose here, the Xi values reside close to xp, then the conclusions 
in the following should hold when the weight factor is applied to the dipole cross section 
depending on ar_i_. A steep increasing as — comes from the behavior of the function 
Ki{r]r±) ~ l/irirx) at this limit. Concerning the dependence on M, as the invariant mass 
increases the contribution looks smaller. 

Regarding the longitudinal contribution, the lower plot in Fig. (0), the weight function 
selects smaller dipole sizes (and 7*g transverse sizes) in comparison with the transverse con- 
tribution. Moreover, the function is narrower as M increases, meaning that larger invariant 
mass scans smaller r±. A well known fact is that the longitudinal contribution is higher 



twist, i.e. it is suppressed by a power of 1/M when compared with the transverse one |]T0 



This feature actually remains in the Drell-Yan case. Moreover, the peaks appearing in the 
plot are due to the balancing between the asymptotic behavior at ^ of the function 
KQ{r]r±) ~ — log(r7r_|_) and the linear r± factor in Eq. ([l^) . 
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0.1 0.2 ^^^^ 0.3 0.4 0.5 

FIG. 2. The longitudinal and transverse contributions for W{r_i_, M'^) as a function of the 7*^^ 
transverse size r± at fixed lepton pair mass M, ior xi ^ xp = 0.525. 

Having addressed the main features of the color dipole framework, namely kinematic 
definitions and the description of Drell-Yan process in the rest frame, in the next section 
one introduces our model for the dipole cross section satisfying unitarity requirements. 

III. THE GLAUBER-MUELLER APPROACH 

The cross section for a color dipole-nucleon scattering is a well known quantity, which was 
first proposed in the BFKL framework [^]. The dipole interacts with the target through a 
perturbative Pomeron, described in terms of the ladder diagrams. From the /cj^-factorization 
framework [^], the scattering process can be written as the convolution of the projectile 
impact factor and the unintegrated gluon structure function from the target, whose dynamics 
is determined by the evolution kernel. The possible orderings in the transverse momentum 
fc-r in these graphs produce the DGLAP or the BFKL dynamical evolutions. In particular, 
considering small r± configurations from the dipole and the fc-r-factorization one obtains. 
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M^,rx) = ^rixG^^^-^^(x,A) (15) 
6 r_|_ 

where xG^'^^^^{x, Q^) is the standard DGLAP gluon distribution at momentum fraction 
X and virtuahty scale = A/r\. An extensive phenomenology has been made using the 
result above for the inclusive structure function and the vector meson production |T^. In 
particular, we call attention to the specific value of the scale Tq appearing in the virtuality 
scale = Tq/t^. We use the Tq = 4 throughout this paper, however other values are 



equivalent at leading logarithmic level [19 



A well defined feature from the data on F2 and on the gluon distribution at high energies, 
i.e. smaller x, is that they present a steep increasing as x decreases. Indeed, experimentally 
F2 ~ xG ~ where the exponent ranges from 0.08 (Regge phenomenology) up to 0.5 (LO 
BFKL calculations). Such a behavior extrapolated to asymptotic energies violates unitarity 
requirements and a control should be considered. The scale where these effects start to 
be important is associated to a region between hard and soft dynamics |^ (pQCD versus 
Regge) or belonging to the high density QCD domain @]. Here, we are interested in the last 
case (for a recent review, see [^). In QCD, the taming of the gluon distribution at high 
energies is taken into account through multiple Pomeron scattering encoded in the eikonal- 



like frameworks [Q. Such a procedure provides the unitarization of the Born Pomeron 
cross section leading to a softer growth with energy. Indeed, the asymptotic calculations 
have produced an unified ln(l/a;) pattern for the cross section and gluon function instead of 
a truly saturated one p5[] . 

In this work we use Eq. (|T5|) , where the standard DGLAP gluon distribution is replaced 
by the modified one. In the following we shortly review how the unitarity corrections are im- 
plemented through the Glauber-Mueller approach [Q. The starting point for the derivation 
is the interaction of a virtual probe particle, in our case a virtual gluon, with the nucleon. In 
the space-time picture of this process, the virtual probe decays in a gluon-gluon {GO) pair 
having transverse separation r^. In the high energy limit, is considered frozen during 
the interaction for x « l/(2m7v-RAr), where the nucleon has mass ttln and geometric 
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transverse size Rn- 

The absorption cross section of a virtual gluon {G*) with virtuahty and Bjorken x 
can be written in the form, 

(x, Q')= f'dzf^ \^^* (Q^ rx, X, z)\' a^ix, r^) , (16) 

Jo J TT J 71 

where z is the fraction of the energy carried by the gluon, b is the impact parameter variable 
and \&'^* is the wavefunction for the transversally polarized gluon generating the pair. The 
cross section of the interaction GG pair with the nucleon depends on energy x and 
transverse separation r^. This description is valid in leading ln(l/x) approximation, however 
in the double log approximation (DLA) of perturbative QCD one obtains 0], 

a^x, rx) = ^rlxG^«^^^(x, ±) . (17) 

The unitarity constraint to the cross section above is expressed by the eikonal-like 
Glauber (Mueller) formula, hence the the gluon structure function can be written as p3|. 



The explicit integration limits for the z (rewritten through the variable x') and transverse 
separation come from the physical kinematic range allowed in the process (for detailed 
discussions, see 0]). The impact parameter b dependence, is parametrized in the profile 
function S{b). It contains information about the angular distribution of the scattering in 
the nucleon case and the nucleon distribution inside the nucleus in the nuclear case. 

The Born term, in the expansion of Eq. (ITM) respect to cr^^, provides the DGLAP 



'N ' 

evolution in double logarithmic approximation (DLA). The remaining terms in the series 
contribute to the saturation effects to the Born term. For simple calculations, the profile 
function S{b) is parametrized as a Gaussian distribution, S{b) = ^;r^e "a , where Ra is the 
target size, which is a free parameter to be determined from data. Then, putting all together 
and performing the integration over impact parameter b in Eq. (|18D, one obtains 

xG(x,Q2) = ^/'^/'^''°^ (^j, + \n[Ka{x',rl)] + E,[Ka{x',rl)]) , (19) 
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where 7^; and Ei{x) are the Euler constant and the exponential integral, respectively. The 
packing factor = {3TTasr'j_/2R\) xG^'^^^^ , sets the scale where saturation effects are 
starting. Namely, the saturation scale Ql is defined through the expression hg{x,QI) = 1. 

Since the Glauber-Mueller approach is valid in DLA, for practical reasons in Refs. |^,|^ 
a procedure was introduced to extend the formalism to the full experimental kinematic 
range available. The final result contains the full DGLAP kernel corrected by contributions 
calculated in DLA, 

xG%'^{x, Q^)= xG{x, Q2) [Eq.{l8)] 

+a;G^^^^^(a;, Q') " ^ £ ^ /J' ^x' G^^^^^ {x' , Q') , (20) 

where this modification is necessary to obtain a realistic approach in the region of not very 
small-x. The above equation includes xG^'^^^^{x, Q^) as the initial condition for the gluon 
distribution and gives xG^'^^'^^{x, Q^) as the first term of the expansion with respect to 
K,G- One needs to subtract the Born term of Eq. (|T^) in order to avoid double counting, 
which is the meaning the last term in equation above. 

From the discussions and definitions above, we should use as dipole cross section in our 
further calculations the following expression, 

afl\x,r^) = ^r'xG%'\x,^). (21) 

The resulting corrected gluon distribution (Eq. (^)) has been applied for a comprehen- 
sive phenomenology in DIS process, considering the formulae above as the gluon input for 
the observables calculated in the Breit frame (structure functions, F2 slope, etc.). Recently, 
the CM dipole cross section has been applied in calculations of the DIS structure functions 
in the dipole color picture (see, Ref. [0). For instance, the structure functions description, 
in particular the latest F2 data, can be seen at figures (3) and (5) of the Ref. p6[, using that 
formalism. 
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FIG. 3. The color dipole cross section as a function of the dipole size r = ar± at fixed X2 = W~^. 
The solid line corresponds to the result without unitarity corrections; the dotted line is the GM 
prediction using 

i?2 ^ iQ and the dashed one is for = 5 GeV ^. 

Now, one discusses in a detailed way the main characteristics emerging from the Glauber- 
Mueller dipole cross section. To do this, in Fig. (^) one shows the Glauber-Mueller dipole 
cross section as a function of dipole transverse size r = ar±_ at fixed momentum fraction 
X2- For sake of a better illustration on the partonic saturation effects, one takes a very 
small value for x = 10~^. Hereafter, one is using the GRV gluon distribution at leading 



order [^. Here we use the GRV94 parametrization, since it has been considered as a 
robust in a comprehensive phenomenology concerning unitarity corrections P,|5|, p6| , |28| , |30| , 
which includes a related R determination, intrinsic to the Glauber-Mueller approach |^ 



The use of the others pdf's |T6|, implies in a determination of the corresponding value of 
the parameter R as well as the enhancement of the already present uncertainty about the 
nonperturbative contribution. For sake of illustration, in order to test the sensitivity to 



the choice of pdf's set [|I6|, in the Ref. fl^ such a study has been performed and it has 
been found that the deviations among the parametrizations are more important in the high 
virtuality region (very small r). However, in this region the color transparency behavior 
~ dominates, absorbing possible sensibility in the specific selected parametrization. 

The solid line corresponds to the standard DGLAP calculation (without saturation), Eq. 
(p!5|), whereas the remaining ones result from unitarity corrections for two different target 
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sizes. The general shape in terms of the dipole size comes from the balancing between the 
color transparency a ~ behavior and the gluon distribution shape. These features are 
depicted in the plots in Fig. where one verifies a visible scaling of xG{x, Q"^) versus 
(left plot) and its dependence on r = ar±_ (right plot). 

The difference in the strength of the unitarity corrections associated with the target size 
is a well known fact. From the data analysis, its value ranges from i?^ = 5 — 10 GeV^^, 
where smaller radius produces strongest corrections [ see Fig. (|)]. For our calculations 
we choose the low value = 5 GeV~^, corroborated by studies in the inclusive structure 



function and its derivative [28 
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FIG. 4. The plot on the left shows the GM gluon distribution (GRV parametrization input) as 
a function of the scale Q"^ = at fixed x = 10^^. On the right, GM gluon distribution versus 
r and the color transparency behavior adip ^ C r"^ (for illustration one uses a free normalization 
C = 15). 



For sake of comparison, one considers the phenomenological model of Ref. (GBW), 
which has produced a good description of HERA data in both inclusive and diffractive 
processes. It is constructed interpolating the color transparency behavior (Jdip ~ t5_ for small 
dipole sizes and a fiat (saturated) behavior for large dipole sizes cr^ip ~ ctq (confinement). 
The expression has the eikonal-like form. 



1 — exp 



(22) 



where Qg = 1 GeV^ and the three fitted parameters are do = 23.03 mb, xq = 3.04 10 
and A = 0.288 and Ro{x) = {x/xqY^"^ is the saturation radius. In GBW, saturation is 
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characterized by the x-dependent saturation radius Q'^{x) = 1/Rq{x) instead of the scale 
coming from Glauber-Mueller, k,g{x,QI) = 1, which can be easily extended for the nuclear 
case [g. 

Although GBW and CM are distinct approaches, the small x DIS data are equally well 
described by both models. In particular the structure functions have been systematically 
described using the GM formalism, see for instance Refs. [|,^,^,^]. The main advanta- 
geous feature of GM in relation to the GBW is the dipole cross section providing a deep 
connection with the gluon distribution, the leading quantity at high energies. Concerning 
GBW, we point out the following shortcomings and disadvantages in comparison with GM 
approach: (a) it is strictly a parametrization available for the small x HERA data; (6) there 
is no direct connection with the gluon content; (c) it does not match DGLAP evolution 
equations; (d) it does not consider the impact parameter dependence of the process; (e) it 
leads to a quite strong saturation scenario in contrast with the other available approaches; 
(/) concerning the hadron-hadron collisions, using GBW to calculate the pion-proton total 
cross section (convoluting the dipole cross section with the pion wavefunction) it predicts 
non-realistic results, i.e. the cross section saturates at ~ 23.03 mb in high energies. 

In Fig. (j^) one shows the Glauber-Mueller dipole cross section as a function of the 
dipole size r = ar± for two typical X2 values. In the lower plot, for X2 = 10~^, the GM cross 
section underestimates the GBW one. However, as X2 decreases the gluon distribution in 
the proton rises and the dipole cross section increases. This feature is depicted in the upper 
plot, for a small X2 = 10~^, where GM overestimates GBW by a factor two at intermediate 
r ~ 0.3. An immediate consequence from the plots is that our prediction lies lower than 
GBW at X2 ~ 10^^ and higher for smaller X2- We discuss these features in a theoretical and 
phenomenological point of view when performing the comparison with available data in the 
next section. 
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FIG. 5. The GM dipole cross section as a function of the dipole size r = arj_ at two typical X2 
values. The GBW result is also shown for sake of comparison. 

IV. RESULTS AND DISCUSSIONS 

This section is devoted to perform theoretical predictions for the available data on 
DY process and the forthcoming ones from RHIC or LHC In the previous section, we 
presented a parameter-free Glauber-Mueller dipole cross section which matches leading log 
gluon evolution and contains corrections from unitarity effects (parton saturation) at higher 
energies. Therefore, this provides a clear dynamical description of the observables depending 
on the gluonic content of the target, also when it is a nuclear one. 

Although perturbative QCD provides reliable results at small distances (small dipole 
sizes), the nonperturbative sector is still far from being completely understood. The usual 
pdf's are evolved from a perturbative initial scale Qq = Mq ^ 1 GeV^, and there is little 
information about the behavior at < Ql, where the perturbative description is not even 
justifiable. In general one makes use of Regge phenomenology to estimate those contributions 
(see, for instance |]19[); and extrapolating to lower virtuality regions (large dipole sizes) one 
needs an ansatz regarding the nonperturbative sector. 

The use of the GRV94 parametrization [^] in our calculations, bearing in mind that = 
meaning its evolution initial scale is = 0.4 GeV^, allows to scan dipole sizes up to 
''cut = ^ GeV~^ (= 0.62 fm). The cut off r(cut) defines the transverse distance scale matching 

18 



the perturbative and nonperturbative sector. For the most recent parametrizations, where 
Qq ~ 1 GeV^ (rcut ~ 0.4 fm) the amount of nonperturbative contribution in the calculations 
should increase. An additional advantage is that GRV94 does not include non-linear effects 
to the DGLAP evolution since the parametrization was obtained from rather large x values. 
This feature ensures that the parametrization does not include sensible unitarity corrections 
(perturbative shadowing) in the initial scale. 

Now, we should introduce an ansatz for the large transverse separation region. A more 
phenomenological way is to match the pQCD dipole cross section with the typical hadronic 
one cr^AT at rcut, for instance as performed in |19|. Nevertheless, due to the large growth of the 
pQCD dipole cross section at high energies and to take a more simple technical procedure we 
choose an alternative way: the gluon distribution is frozen at scale rcut, namely xG{x, Qcut)- 
Then, the large distance contribution r > rcut reads as, 

aS^(x, r > rcut) = ^ dt ^ ^) • (23) 

"J '"cut 

In a more rigorous analyses, one should substitute the freezening scale Qcut f*^^ ^^e saturation 
scale Ql{x) to take into account a realistic value of the gluon anomalous dimension in all 
kinematic regions. 

Concerning our ansatz for the large r region, one verifies in Ref. that it produces 
reasonable data description, mainly the normalization of the structure functions. An im- 
provement for the non-perturbative contribution is performed there, where a cut off in the r 
integration (0 < r < Vcut) and the addition of a soft Pomeron term are considered. However, 
if we do not introduce this improvement, the normalization of the structure functions re- 
mains unaffected. For completeness, the consideration of a soft term for the non-perturbative 



contribution was also taken into account in Ref. |]T9[. However, for the same reasons above, 
we choose the technically more simple procedure. 
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FIG. 6. The profile HT,L{r±, M'^) as a function of the 7*5 transverse separation r± at typical 
mass M values. One uses X2 = and GRV parametrization input. 

To illustrate the role played by the small and large transverse separations in the de- 
scription of the observables, in Fig. (^) one shows the profile of the r^-integration from 
Eq. ( ]I3| ) as a function of the 7*g transverse separation r± at typical mass M values. It is 
labeled here as HTL{r±,M'^). The momentum transfer is fixed at X2 = 10~^, since the low 
X2 data available lie at this magnitude. For the proton structure function M^), 
describing the quark content of the projectile, we use the ALLM updated parametrization 
||3l| (good agreement with HERA data at large x) . Both transverse and longitudinal profiles 
are presented. 

The main contribution for the profiles comes from the asymmetric peaks which are shifted 
to larger r± as M diminishes. For instance, in the transverse profile, the peak lies at 
r± ^ 0.06 fm for M = 10 GeV whereas at M = 5 GeV it takes values r± ^ 0.1 fm. As 
it was verified in Sec. (2), the longitudinal sector selects smaller transverse sizes r± than 
the transverse one. Indeed, from the upper plot, non-zero contributions are obtained up to 
large r± ^ 1 fm. The higher twist character of the longitudinal piece is verified through the 
magnitude scale of HL{r±, M^). Due to the fact that the nonperturbative sector dominates 
for Tcut = {c(r±)cnt > 0.62 fm using the GRV input, a significant part of the contributions 
comes from the soft region where a is small (soft quark). This is in agreement with the 
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expectations that important soft contributions take place in Drell-Yan process (see related 



discussions at ). 

Now, we are able to compare the results with the available data. Since the color dipole 
picture is valid at small momentum fraction X2, one needs to select the experimental data 
covering this requirement. The lowest X2 data were obtained in the fixed-target dimuon 



production at E772 Collaboration |[18|, where we select the points with X2 < 0.1 following 



the similar procedure of [0. In Fig. (|^, one presents the calculation Eq. ([13|) using the 



Glauber-Mueller dipole cross section (the solid line) at fixed xp and center of mass energy 
^/s = 38.8 GeV (0.03 < X2 < 0.09). It should be stressed that this kinematical region scans 
the validity limit of the color dipole picture. The curves underestimate similar calculations 
in Ref. Il^], which uses the phenomenological GBW dipole cross section. Such a result is 
actually expected from our conclusions in the previous section where GM underestimate 
GBW at X2 = 10~^ [see lower plot in Fig. (H)]. 
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FIG. 7. The DY differential cross section per nucleon versus M for the available small X2 data 
[18] {^/s = 38.8 GeV) at fixed Xi^' in pD reaction. The solid line corresponds to the Glauber- Mueller 
dipole cross section. The dashed one is the LO Breit system calculation and the dot-dashed line 
corresponds to the Glauber-Mueller added of a reggeon contribution. 

The experimental data analysed above are for x > 10~^. Therefore, the Eq. (|2lD is 
known to have corrections at larger x values due to the exchange of quarks rather than 
gluons, in the t-channel, corresponding to a reggeon instead of a Pomeron exchange ||32|| . 
The secondary reggeon contribution corresponds to an amplitude with quark-antiquark pair 
t-channel exchange. The leading double-logarithmic asymptotics of such an amplitude was 
calculated in perturbative QCD in Ref. |Q. The quark-antiquark cut occurs in the j-plane 
at ujoit) = ^2Cpas/iT, where Cp = {N^ — l)/2Nc and is the strong coupling constant. 
That value is very close to the phenomenological intercepts of the u, p trajectories, i.e. 
aR^O) ~ 0.5. Our expression for the dipole cross section, Eq. (^if) , considers only sea 
contribution for the process (gluon radiation), being equivalent in the Regge terminology 
to the hard Pomeron. The correspondent valence-like, which corresponds to the reggeon 
contribution, is lacking in our analysis above. In order to simulate the valence content in 



the calculations, we parametrize that piece in the following form [34]: 

a,n{x,r)=Njnr'x'^-''''{l-xf, (24) 
where we have used the reggeon intercept ^/^(O) = 0.5475 and the threshold factor for the 



large x region ^J]. To reproduce similar results as presented in Ref. [0, one considers 
the constant value Njr = 8 (to obtain a ct/r in GeV~^). The factor ensures the correct 
scaling. 

In the plots (Fig. |^ one shows also the LO Breit system calculation, Eq. (|I]), which 
is the dashed line. The color dipole result considering only the gluon content (sea quarks), 
Eq. (pTD, lies below LO fast proton frame one at X2 ~ 10^^ (where the presented data are 
available). However, in this kinematical region the valence quark content competes with the 
sea one and such a difference should be expected. Now, we introduce the valence content, 
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parametrized in Eq. (^4]) and added to the Eq. (pT]): as a result our data description has 



been improved, equivalently to the calculations in considering the GBW model. As 
X2 decreases the gluonic content of the target drives the observables and the color dipole 
considering only Eq. (|2T| ) should produce quite reliable results. We have verified this feature 
and have found that the reggeon contribution for the RHIC energies is completely negligible. 

Here, some comments about higher order corrections are timely. The color dipole ap- 
proach results for total (virtual) photon cross section are equivalent to those ones obtained 
by /cy-factorization in the leading logarithmic approximation. However, the inclusion 



of higher order effects in the fcr-factorization approach turns the equivalence incompatible: 
the conservation of the transverse positions and sizes of the colliding objects is violated . 



Therefore, the introduction of higher orders contributions into the dipole cross section must 
be taken with some care. Moreover, deep inelastic and Drell-Yan have a quite different sce- 
nario concerning NLO and NNLO corrections. In DIS, calculations considering up to NNLO 



resummations have been performed and it was found that they are small |3^. Instead, in 
Drell-Yan even the NLO calculations produce corrections up to a factor of two, diminishing 
as the energy increases |0,[TB|. Keeping in mind the discussion above, at the moment we are 
unable to perform in the dipole color picture an equivalent NLO (Breit frame) calculation, 
since at the present the wavef unctions are not still available at NLO accuracy. 

In order to address the color dipole picture at high energies, the DY differential cross 
section for RHIC energies, a/s = 500 GeV, is shown in Fig. (||) for the same fixed xp- There, 
the X2 reaches values of order 10~^ and unitarity effects are important. The solid lines are the 
Glauber-Mueller estimates, Eq. (0), and the dot-dashed ones are the rest frame calculations 
with DGLAP gluon distribution, Eq. ([T5|). The curves overestimate similar calculations in 
Ref. |T^, which uses the phenomenological GBW dipole cross section. Such a behavior is 
expected from our previous conclusions where GM overestimates GBW at smaller X2 due 
to the growth of the gluon distribution at higher energies [see the upper plot in Fig. (|^)]. 
Concerning the rest frame non-corrected DGLAP input, the Glauber- Mueller underestimates 
them due to the significant corrections coming from unitarity effects (parton saturation). 
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Moreover, we have obtained a result almost similar to the LO Breit frame calculations at 
the RHIC energies, suggesting a good consistency in both frameworks. From the plots one 
verifies that the deviations are more significant as M diminishes, corresponding to smaller 
X2 values. In absolute values, the corrections at RHIC energies reach up to ~ 20 %. 
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FIG. 8. The DY differential cross section per nucleon versus M for the RHIC energies {y/s = 500 
GeV) at fixed xp in reaction. The solid line corresponds to the Glauber- Mueller dipole cross 
section whereas the dot-dashed one is the non-corrected DGLAP calculation. 



As final comments, we address additional advantages of the color dipole picture in the 
DY case. For example, it allows to obtain the transverse momenta px distribution for the 



process already at leading order calculation [0. Instead, in the parton model the lepton 
pair has no transverse momentum due to the assumption that in the partonic subprocess the 
longitudinal momenta are bigger than the transverse ones (partons are collinear). Therefore, 
an alternative way to solve this trouble is to introduce an intrinsic pt for the initial state 
interacting partons. However, such an assumption is not sufficient to describe the measured 
Pt distributions. Considering the Compton and the annihilation subprocess the leptons 
acquire transverse momentum and the p-p-dependence can be calculated in pQCD. The 
resummations produce ln^(m^/p^) terms which are large as the transverse momentum 
goes to soft values — ^ (the perturbative expansion breaks down). Thus, the color 
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dipole description is a nice tool to calculate that distributions since the above difficulties are 
avoided. One intends to address carefully this issue in a next calculation. 



V. CONCLUSIONS 



The Drell-Yan is an important process testing the quark (antiquark) content of the 
hadron target. The measured observables are Lorentz invariant, whereas the parton descrip- 
tion is frame dependent. Calculations in the fast proton system have provided a perturbative 
understanding of DY up to higher orders. On the other hand, the color dipole picture allows 
a simple description of DY driven by the gluonic content (sea quarks) of the target. The 
quark from the projectile radiates a photon decaying into a lepton pair. The basic blocks in 
the color dipole are the LC wavefunctions and the dipole cross section. The former is cal- 
culated from perturbation theory and the latter one is modeled taking into account general 
properties of both hard and soft domains. 

We have found that the LC wavefunction for the 7*g configuration plays the role of 
a weight function for the different transverse separations r± (as well as dipole sizes ar±) 
contribution to the process. Namely, small transverse separations are selected by both the 
transverse and longitudinal pieces. However, the transverse contribution can select non- 
negligible large sizes r±. In addition, larger invariant mass M scans smaller ■j*q separations. 
Moreover, the longitudinal piece is higher twist (suppressed by a factor 1/M^). These 
features are also present in the deep inelastic case due to the similarity between the LC 
wavefunctions expressions. 

Concerning the dipole cross section, here we consider the Glauber-Mueller approach, 
which takes into account the corrections to the standard DGLAP formahsm due to the 
unitarity requirements. The taming of the parton distributions (parton saturation) at high 
energies is performed considering the multi-scattering assumption from the Glauber-like 
(eikonal) formalism. We have found a distinct behavior at both low and large X2 when 
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performing a comparison with the phenomenological GBW modeL The main source of the 
deviation is that GM depends on the gluon distribution, which increases as X2 diminishes. 
These features produce distinct results at current energies and in the forthcoming mea- 
surements. An important verification is that a non-neghgible amount of non-perturbative 
contribution is present in the cross section. Although the LC wavefunctions suppress large 
transverse separations, a large cross section at small X2 compensates the suppression pro- 
ducing significant soft content. 

The current low X2 available data lie in values ranging from 0.03 < X2 < 0.09, actually 
testing the validity limit of the color dipole picture. Our results considering only GM dipole 
cross section underestimate the experimental measurements since color dipole includes only 
the sea quark content (gluon radiation) from the target. In the realistic case, for this 
kinematical region the valence and sea quarks have both a significant contribution in the 
cross section. We have parametrized the valence content through a reggeon exchange and 
the results turn out equivalent to the data description claimed in Ref. 11^. It was found 
also that issues related to higher orders contributions in the color dipole picture should be 
taken carefully. 

As the energy reached in the forthcoming experiments increases, the saturation effects 
should turn out to be more relevant. We perform estimates for the RHIC energies and have 
found that the unitarity corrections are important in the description of the cross section. 
We expect that such correction should be larger at LHC, since the X2 values probed there 
would be smaller than in RHIC. 

The quite simple scenario for DY process in the rest frame allows to extend the approach 
to the nuclear case and also get information on pt distribution. The higher energies soon 
available will demand a well established knowledge on the nuclear gluon distribution which 
can be the input for the nuclear dipole cross section in the color dipole framework. This 
approach should be a useful tool to perform pQCD estimates for the future experimental 
measurements. 
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